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Abstract: A common quark potential that captures the essential traits of the QCD quark-gluon dynamics is 
expected to (i) interpolate between a Coulomb-like potential (associated with one-gluon exchange) and the 
infinite wall potential (associated with trapped but asymptotically free quarks), (ii) reproduce in the intermediary 
region the linear confinement potential (associated with multi-gluon self-interactions) as established by lattice 
QCD calculations of hadron properties. We first show that the exactly soluble trigonometric Rosen-Morse 
potential possesses all these properties. Next we observe that this potential, once interpreted as angular 
momentum dependent, acquires a dynamical 0(4) symmetry and reproduces exactly quantum numbers and level 
splittings of the non-strange baryon spectra in the SU (2)/ £g) 0(4) classification scheme according to which 
baryons cling on to multi-spin parity clusters of the type ( -j, y) <8 [(5, 0) ® (0, |)], whose relativistic image is 
VV- ma • Finally, we bring exact energies and wave functions of the levels within the above potential and thus put 
it on equal algebraic footing with such common potentials of wide spread as are the harmonic-oscillator- and the 
Coulomb potentials. 



PACS: 02.30.Gp, 03.65.Ge,12.60.Jv. 

1 Introduction 

The non-strange baryon spectra below ~ 2500 MeV reveal, isospin by isospin, as a striking phenomenon mass 
degenerate series of K pairs of resonances of opposite spatial parities and spins ranging from ^ to (K —5) 
which terminate by a highest spin-(if + |) resonance that remains unpaired 1 1 1. Such series (displayed in Fig. 1) 
perfectly fit into SU(2)i ® 0(4) representations of the type If) ® i{h> ^) ® §)]' an observation due again 
to Q). The appeal of the above classification is twofold. On the one side, up to the A(1600) state which is likely to 
be a hybrid, no other resonances drop out of the proposed scheme. Also the prediction of less "missing" resonances 
relative to others schemes falls under this issue. On the other side, due to the local 0(4) ~ 0(1, 3) isomorphism, 
the non-relativistic 0(4) multiplets have as an exact relativistic image the co variant high-spin degrees of freedom 
given by the totally symmetric rank-J\" Lorentz tensors with Dirac spinor components, ijj f _ ll ... f _ lK known as spin- 
{K + 5) Rarita-Schwinger fields |2|. In this fashion, one can view the series of mass degenerate resonances of 

alternating parities and spins rising from | to [K + as rest frame t/Vi-. ^k- °f mass m an d l°°k for possi- 
bilities to generate such multiplets as bound states within an appropriate quark potential. Although the degeneracy 
of the non-strange baryons follows same patterns as the states of an electron with spin in the Hydrogen atom, the 
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Figure 1: Experimentally observed baryon resonances (l.h.s.) N and (r.h.s.) A. The dash-point lines represent 
the series mass average. Notice that the resonances with masses above 2000 MeV are of significantly lower 
confidence but those with masses below 2000 MeV where the degeneracy is very well pronounced. Empty squares 
denote predicted ("missing") states. Typical, the systematical lack of a parity partner to the first highest spins 

Djs, Fir, an d Hi ii (the last two being among the "missing" N states). 



level splittings are quite different. The mass formula that fits the N(A) spectrum has been encountered in Ref. ||3) 
as 
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V ~ 4 ' 
/a = 600 MeV, g N = 70 MeV, 



a = K + l, I = N,A, 



g A = 40 MeV . 



(1) 

(2) 



and contains besides the Balmer-like term, (~ — 1/<t 2 ), also its inverse. In effect, the baryon mass splittings 
increase with a. The degeneracy patterns and the mass formula have found explanation in Ref. 1 4 1 within a version 
of the interacting boson model (IBM) for baryons. To be specific, to the extent QCD prescribes baryons to be 
constituted of three quarks in a color singlet state, one can exploit for the description of baryonic systems algebraic 
models developed for the purposes of triatomic molecules, a path pursued by Refs. |5 1. An interesting dynamical 
limit of the three quark system is the one where two of the quarks are "squeezed" to one independent entity, a di- 
quark (qq), while the third quark (q) remains spectator. In this limit, which corresponds to U (7) — ► U (3) x U(4), 
one can exploit the following chain of reducing U(4) down to 0(2) 



17(4) D 0(4) D 0(3) D 0(2) , 

N K I m 

K = N,N- 2, ...1(0), l = K,K- 



1. 



,0, 



< I. 



(3) 
(4) 



in order to describe the rotational and vibr ational (rovibron) modes of the (qq) — q dumbbell. In so doing, one 
reproduces the quantum numbers describing the degeneracies in the light quark baryon spectra by means of the 
following Hamiltonian: 



, K K 
C 2 (so(4)) ( -. - 



H-H a = -//(4C 2 (so(4)) + 1)" 1 
(K + l) 2 - 1 (K K 



giC 2 (so(4)). 



(5) 
(6) 
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with C2(so(4)) being the second so(4) Casimir operator. In the second row of Eq. we indicate the quantum 
numbers associated with the respective group of the chain. Here, N stands for the principle quantum number of 
the four dimensional harmonic oscillator associated with U(4), K refers to the 0(4) four dimensional angular 
momentum, while I, and m are in turn ordinary three- and two angular momenta. In Ref. |4| the interested reader 
can find all the details of the algebraic description of the nucleon and A resonances within the rovibron limit. 

Yet, as a principle challenge still remains finding a suitable quark potential that leads to the above scenario. In 
the present work we make the case that the exactly soluble trigonometric Rosen-Morse potential is precisely the 
potential we are looking for. 

The paper is organized as follows. In the next section we analyze the shape of the trigonometric Rosen-Morse 
potential. In section III we present the exact real orthogonal polynomial solutions of the corresponding Schrodinger 
equation. The paper ends with a brief concluding section. 



2 The shape of the trigonometric Rosen-Morse potential 

We adopt the following form of the trigonometric Rosen-Morse potential [6 1,|7 1 

v(z) = -2b cot z + a(a + 1) esc 2 z , a > -1/2, (7) 

displayed in Fig. 2. Here, 

z= r -, v(z) = V(z)/(H 2 /2md 2 ), e n = E n /(h 2 /2md 2 ) , (8) 

the one-dimensional variable is r = |r|, d is a properly chosen length scale, V (r) is the potential in ordinary coor- 
dinate space, and E n is the energy of the levels. Our point here is that v(z) interpolates between a Coulomb-like 
and an infinite-wall potential going through an intermediary region of linear, and quadratic (harmonic-oscillator) 
dependences in z. To see this (besides inspection of Fig. 2) it is quite instructive to expand the potential in a Taylor 
series which for appropriately small z, takes the form of a Coulomb-like potential with a centrifugal-barrier like 
term (if a were to be a positive integer) provided by the esc 2 z part, 

, . 2b a(a + 1) 

v(z) « + 1 2 ' • (9) 

z z A 

In an intermediary range where inverse powers of z may be neglected, one finds the linear plus a harmonic- 
oscillator potentials 

, , 2b a(a + l) , 
v(z)tt— z + \ > z 2 . (10) 

A 15 

Finally, as long as cot z — ► 00 and esc z — > 00, the potential obviously evolves to an infinite wall. The above 
shape captures surprisingly well the essentials of the QCD quark-gluon dynamics where the one gluon exchange 
gives rise to an effective Coulomb-like potential, while the self-gluon interactions produce a linear confinement 
potential as established by lattice calculations of hadron properties. Finally, the infinite wall piece of the trigono- 
metric Rosen-Morse potential provides the regime suited for trapped but asymptotically free quarks. By the above 
considerations one is tempted to conclude that the potential under consideration may be a good candidate for a 
common quark potential of QCD traits. 

In the next section we present the exact solutions of the Schrodinger equation with the trigonometric Rosen-Morse 
potential. 
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Figure 2: The proximity of the (~ cot z)- to the (~ -) potential. 



3 Energies and wave functions of the levels within the trgonometric Rosen- 
Morse potential 



The three-dimensional Schrodinger equation with the trigonometric Rosen-Morse potential (tRMP) reads: 

V 2 ^(z)+ (lb cot z - a( " t ^ + e ) 4>(z) = , 
\ sin z I 



(11) 



and is solved in polar coordinates in the standard way by separation of variables. In effect, the wave function 
factorizes according to 

V ,( z ) = R&Yffrfl, 1 = 0,1,2,..., \fi\<l, (12) 
where Y^iO, <j>) stand for the standard spherical harmonics, and R(z) satisfies the one-dimensional equation 

^ R(z) +(2bcotz- a(a -± 1} + e) R(z) = 



sm z 



(13) 



This equation (up to inessential notational differences) has been solved in our previous work |8|. There, we 
exploited the following factorization ansatz 



R{z) = e- az/2 {l + cot 2 z)^ 1 C^> a) {cot z) , 



(14) 



with a and (3 being constant parameters. Upon introducing the new variable x = cot z, substituting the above 
factorization ansatz into Eq. J14i . and a subsequent division by (1 + 

!B 2)(l-/J)/2 one 

finds 



2 ,d 2 C^-^ix) /a a 
{1 + X) d* +2 \2 +fiX 



dCW- a Hx) 
d x 



{-ail -/3) + 2b)x + ((f) 2 -(l-/3) 2 +e m ) 
+ ((-/3(l-/?)-a(a + l)) + J -)c^ a \x)=Q. (15) 
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This equation is suited for comparison with 

(l + x 2 ) 7% W +2(^+/3x) d/ ^ -m(2/3 + m-l)^' Q) W = 0, (16) 

which being of the form of the text-book hypergeometric equation |9 1,| 1(31. 11 II can be cast into the self-adjoined 
Sturm-Liouville form given by 

, {x f2j^ + j_(^^yj^ + ^ n ^ {:c) . , (17) 

or \ ax J ax 

s{x) = l + x 2 , w^^ a) (x) ^ s{x) f3 - 1 e' acot ' lx , X m = -m(20 + m-l), -oo < at < oo. (18) 

However, while the standard textbooks consider exclusively s(x) functions which are at most second order poly- 
nomials of real wots, in which case 

w(a)s(a)x l = w(b)s(b)x l = , VZ = integer, (19) 

holds valid, the roots of s(x) in Eq. (116b are imaginary. In the former case it is well known that 

• TZ-m {x) would be polynomials of order to, 

• A m would satisfy 

( v dTl{ M (x) 1. i,d 2 s(x)\ 
A ™ = -m\K x i 7 _^ + -( m -l)- 7 ^i , 



g? x 2 d x 2 



(20) 



with K m being the TZm' 01 ^ (x) normalization constant, 
• the first order polynomial would be defined as 



n (M (x) = ^_(ds^pxY (21) 



Kiw(x) \ dx 

• the latter relation would generalize to any m via the so called Rodrigues formula 



1 d m 

^' a) ^ = ^)^^ S ^ m )' (22) 

• w{x) would be the weight-function with respect to which the IZm {%) polynomials would appear orthog- 
onal. 

Within this context the question arises whether the imaginary roots of s(x) = (1 + x 2 ) in Eq. d!6i would prevent 
the TVm (x) functions from being real orthogonal polynomials. The answer to this question is negative. It can be 
shown that also in the latter case 

• the TZm (x)'s are polynomials of order m, 

• they can also be constructed systematically from a Rodrigues formula in terms of the respecified weight 
function, 

• but only a finite number them will be orthogonal due to the violation of Eq. ( 1191 . 
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From the historical perspective, Eq. Jl 8I > has first been brought to attention by the celebrated English mathematician 
Sir Edward John Routh in Ref. 1 12 1 (modulo the unessential difference in the argument of the exponential from 
the present cot -1 to Routh's tan -1 ), the teacher of J. J. Thomson and J. Larmor, among others famous physicists. 
Routh observed that the weight-function of the Jacobi polynomials, Pm' v \x), takes the form of Eq. (1 1 8i upon 
the particular complexification of the argument and the parameters according to /i — > 77 = a + ib,v — > rf , 
and x — * ix. From that he concluded that Pm' v (ix) is a real polynomial (up to a global phase factor). Later 
on, in 1929, the prominent Russian mathematician Vsevolod Ivanovich Romanovski, one of the founders of the 
Tashkent University, studied few more of their properties in 1 13 1 and it was him who observed that only a finite 
number of them appear orthogonal. While the mathematics literature is familiar with such polynomials 11141 . Ill 51 . 
Ill 61 . 1171 where they are referred to as finite Romanovski polynomials 1 18 1, or, Romanovski-Pseudo- Jacobi poly- 
nomials 1191 . a curious omission from all the standard textbooks on mathematical physics |9|,|10| takes place. 
This might be related to circumstance that the physics problems which call for such polynomials are relatively 
few. Recently, it has been reported in the peer physics literature |8 1, |20| that the Schrodinger equation with the 
respective hyperbolic Scarf and trigonometric Rosen-Morse potentials is resolved exactly precisely in terms of the 
Romanovski polynomials. Moreover, the latter are also relevant in calculation of gap probabilities in finite Cauchy 
random matrix ensembles 1121 1 . In the following, we shall adopt the notion of Routh-Romanovski polynomials for 
obvious reasons. 



3.1 The exact spectrum 

Back to Eq. H5\ we observe that if it is to coincide with Eq. dl6> then the coefficient in front of the l/(x 2 + 1) 
term has to nullify. This restricts the C^ ,a '(x) parameters in the Schrodinger wave function to be 

-J - (l-/3) 2 + e = 0. (23) 
With that Eq. dl5> to which one has reduced the original Schrodinger equation amounts to 



(! + +2 (f + *) + (-,(1 - ffl - a(a + l)) C <-»(,) = 



(24) 



The final step is to identify the constant term in the latter equation with the one in Eq. d!6i which amounts to a 
third condition 

- 0(1 - 0) - a(a + 1) = -m(20 + m - 1) , (25) 

which introduces the dependence of the C^' a \x) functions on the index to, i.e. C^' a \x) — > Cm'"'(i). 
Remarkably, Eqs. \22>\ and d25l > indeed allow for consistent solutions for a, 0, and e and given by (upon renaming 

to by (n — 1)): 



2b 

n + a 



(n + a) 2 - j- ; , n = m + 1 , (26) 



n = -{n + a) + 1 , 
(n + a) 2 

now with n > 1. In this way Eq. J26l > provides the exact tRM spectrum. In effect, the polynomials that define the 
exact solution of the Schrodinger equation with the trigonometric Rose-Morse potential turn out identical to the 
Routh-Romanovski polynomials however with indices that depend on n. As we shall see below, this circumstance 
will become of crucial importance in allowing for an infinite number of orthogonal polynomials (as required by 
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the infinite depth of the potential) and thus for avoiding the finite orthogonality of the bare Routh-Romanovski 
polynomials. 

With that all the necessary ingredients for the solution of Eq. (I24> have been prepared. In now exploiting the 
Rodrigues representation (when making the n dependence explicit), 

1 H n ~ 1 / \ 

C^\x) = 7e-)(*) = Knwi ^ ){x)dx ^ «(aO n_1 ) , ill) 

allows for the systematic construction of the solutions of Eq. J24I . Notice that in terms of w^ n,an ^ (x) the wave 
function is expressed as 



^(cot- 1 x) = V / w (-(" +a )+ 1 ^) (.x)^ ( r ( " +a)+1 ^ ) (x) . (28) 
Next one can check orthogonality of the physical solutions in z space and obtain it as it should be as 

R n (z) Rn> (z) dz = 5 nn > , (29) 



The orthogonality of the wave functions R n {z) implies in x space orthogonality of the 1Z^ n ' a "' (x) polynomials 
with respect to ui^" ,Q "^(x)^| due to the variable change. As long as dco ^ x = —1/(1 + x 2 ) = —l/s{x) then 
the orthogonality integral for the polynomials takes the form 



dx 
s(x) 



yJwW* ^ {x)K<£ » •«»> (x) yjw^ ) {x)U^' ] (x) = S n n , . (30) 



The existence of an infinite number of orthogonal Routh-Romanovski polynomials was made possible on cost of 
the n dependence of the parameters which emerged while converting the Schrodinger equation to the polynomial 
one. 



3.2 The degeneracy in the spectra 



Inspection of Eq. i26i reveals existence of an intriguing degeneracy in the tRMP spectrum. In order to see it let 
us assume that the a-parameter in Eq. H3i takes only integer non-negative a = 0, 1, 2, ...-values. In such a case, 
one immediately reads off from Eq. i26\ that the energy levels for any a = (m + 1 + a) with a = 1,2, 3, ... are 
(T-fold degenerate as a can take all the values between and (a — 1) according to a = 0, 1, ...(a — 1) (see Fig. 3). 
Comparison of the a-degeneracy to the non-strange baryon spectra in Fig. 1 and Eqs. @ is suggestive of the idea 
to interpret the non-negative integer a-values as angular momenta and view the esc 2 (5) term as a non-standard 
centrifugal barrier 

a(a +1) esc 2 (J) — > -^TTTT' = 1 = 0,1,2,.... (31) 
\dJ sm 2 (^) 

In terms of a the mass formula in Eq. (|2) translates to 

4(M g , / -A/ M + l g/ ) 6 2 4fj 

> ( a = 9 + er , o = . (32) 

91 o 91 

Non-standard centrifugal barriers of various types have been frequently exploited in the calculation of the spectral 
properties of collective nuclei. Specifically, in Ref. 1221 use has been made of an angular-momentum dependent 
potential originally suggested by Ginocchio in Refs. |23|. The non-standard centrifugal barrier in this potential 
asymptotically approaches for certain parameter values the physical centrifugal barrier, 1(1 + l)/r 2 while for 
another set of parameters it becomes the Poschel-Teller potential. In our case, for small arguments the esc 2 term 
also approaches the physical centrifugal term as evident from Eq. l|9} and visualized by Figs. 4. Non-standard 
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Figure 3: Degeneracy of energy levels of same a but different angular momenta in I dependent trigonometric 
Rosen-Morse potentials. The curves correspond to b = 60, a value fitted to the TV spectrum. 



centrifugal barriers have the property to couple various multipole modes in nuclei, an example being given more 
recently in Ref. |24|. From now onward we shall adopt non-negative integer values for the a parameter and view 
the esc 2 term as a non-standard centrifugal barrier according to 



2 



V ^=- 2BCOt (d) + ^ l + V C8C2 {d)- ^ 

In so doing we are defining a new angular momentum dependent potential, V/(r), that possesses the dynamical 
0(4) symmetry. Notice that this does not contradict the statement of Ref. |25| according to which only pure or 
screened Coulomb-like potentials are 0(4) symmetric as the theorem of Ref. [25 1 refers to potentials with the 
standard centrifugal barrier only. The b parameter in Eq. (1321 now relates to the potential parameter B as 

(34) 

Next we shall bring down the a index, suppress the b index and change notations according to 

^ (O=0 ' 6) ©^^©' ° = n + l- (35) 

The single-particle wave functions within the angular dependent trigonometric Rosen-Morse potential are straight- 
forwardly calculated from Eq. (1281 . Below we list the first three levels for illustrative purposes: 

• ground state a = 1: 



Is: R W ( V - 



dl V (1 



b(b 2 + l) _ b( r) . /r\ 



8 




Figure 4: The cot (^) potential with the physical centrifugal barrier (left) and the non-standard one (right). 

• first excited state, 

2s: j 
2p: j 

• second excited st 

3s: tfao©) 

3 P : R 31 Q 
3d: R 32 Q 

In Figs. 5 we display as an illustrative example the wave functions for the first two a levels. 
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Figure 5: Wave functions for a = 1, 1 = (left) and a = 2, 1 = 0, 1 (right). 

4 Discussion and concluding remarks 

In this work we made the case that the trigonometric Rosen-Morse potential with the esc 2 term being reinterpreted 
as a non-standard centrifugal barrier provides quantum numbers and level splittings of same dynamical 0(4) 
patterns as observed within the SU (2)/ x 0(4) classification scheme of baryons in the light quark sector. Due to 
local 0(4) ~ 0(1, 3) isomorphism, the potential (-j, If) ® i(h> ^) ® \)\ l eve ls have as a relativistic image 
the covariant field theoretical high-spin degrees of freedom, t/j f _ ll ... f _ lK ■ 

We presented exact energies and wave functions of a particle within the above potential and in so doing put it on 
equal algebraic footing with the harmonic oscillator and/or the Coulomb potentials of wide spread. In this fashion, 

• an element of covariance was brought into the otherwise non-covariant potential picture, 

• the algebraic Hamiltonian in Eq. (|6j describing the 0(4) degeneracy in the N and A spectra was translated 
into a potential model of same dynamical symmetry. 

The 0(4) degeneracy of the N and A spectra seem to speak in favor of quark-diquark as leading configurations 
of resonance structures. Yet, form factors are known to be more sensitive to configuration mixing effects and may 
require inclusion of genuine three quark configurations. As long as the tRMP shape captures the essential traits 
of the quark-gluon dynamics of QCD, we here consider it as a promising candidate for a realistic common quark 
potential that is worth being employed in the calculations of spectroscopic properties of non-strange resonances. 
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